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A fermion model with random on-site potential defined on a two-dimensional square lattice with 
\ TT-flux is studied. The continuum limit of the model near the zero energy yields Dirac fermions 

. with random potentials specified by four independent coupling constants. The basic symmetry of 

the model is time-reversal invariance. Moreover, it turns out that the model has enhanced (chiral) 
symmetry on several surfaces in the four-dimensional space of the coupling constants. It is shown 
that one of the surfaces with chiral symmetry has Sp(n)xSp(n) symmety whereas others have 
U(2n) symmetry, both of which are broken to Sp(n), and the fluctuation around a saddle point is 
QQ , described, respectively, by Sp(n)2 WZW model and U(2n)/Sp(n) nonlinear sigma model. Based on 

these results, we propose a phase diagram of the model. 



Ch ' PACS: 72.15.Rn, 71.23.-k, ll.lO.-z 



I. INTRODUCTION 



Anderson localization has attracted renewed interest since a plenty of new universality classes were discoverd 
One is a system with particle-hole symmetry studied by Gade jl]. Although generic disorderd systems in 
two dimensions are believed to be insulators [Q, she showed that such a system has a random critical point at the 
band center. Typical examples studied so far are the random flux model [^|-^ and the random-hopping model with 
Q TT-flux l^j^j^J^. This class is often referred to as AIII and BDl in the case with broken and unbroken time reversal 
symmetry, respactively, which corresponds to the chiral Gaussian unitary and chiral Gaussian orthogonal ensamble 
I ' of the random matirx theory in the zero dimensional limit ||^ . 
^ Other examples of new universality classes are found in disordered d-wave superconductors described by a mean- 

^\ field Bogoliubov-deGennes Hamiltonian (ll-23|. Altland and Zirnbauer |^,|| classified them by two basic symmetries. 



i.e, spin-ratation and time-reversal into four kinds of classes. With (without) spin-rotation invariance, disordered 
d-wave superconductors belong to class C and CI {D and -Dili), respectively, in the case with broken and unbroken 
time-reversal symmetry. The beta function of the renormalization group for corresponding nonlinear sigma models 
in two dimensions tells that the classes C and CI have localized states only, whereas the classes D and I?III may 
, have dclocalized states. What is particularly interesting is that some of them allow in the nonlinear sigma model 
description the topological term [p^|jl^ , |2^ which plays a crucial role, for example, in quantum Hall transitions [p4[ , 
or the Wess-Zumino-Witten (WZW) term |17 1^,|2^ which gives rise to a nontrivial strong-coupling fixed point. 



Therefore, disordered d-wave superconductors are expected to have quite rich phase diagram. 
, These developments can shed new light on the localization problems. For example, the question whether the 
^ ' density of states of the random fiux model diverges had been a long-standing problem. Recognizing, however, that 
Q it belongs to the Gade's class of ^III type, Furusaki concluded that the DOS actually diverges. Another example is 
O the disordered fermion model with 7T-fiux we will revisit in this paper. This model was already studied by Fisher and 
J> , Fradkin in 1984. They showed that the effective field theory of the lattice model with random on-site potentials 
' belongs to the universality class of 0(2n, 2n)/0(2n)xO(2n) nonlinear sigma model, concluding thereby that all states 
r\ \ are locafizaed. 

' A point here is that the continuum limit of the model includes random potentials specified by four independent 
coupling constants. Since Fisher and Fradkin studied only the case where all the coupling constants are equal, we can 
expect richer phase diagram in the whole space of the coupling constants. Actually, although the generic points in 
this space should fiow into the localized phase mentioned-above, the model has enhanced symmetry in special cases 
and accordingly renormalization group flows in them are toward different fixed points. Another point in this paper is 
concerned with chiral symmetry. It was shown that generic disordered d-wave superconductors in class CI, described 
by Dirac fermions, yield no WZW term ||l7| , ^ although those Dirac fermions have chiral Sp(n)xSp(n) symmetry. 
The reason may be that chirality is not unique in the d-wave Dirac fermions, since the lattice fermion gives the species 
doubling in the continuum limit. Here, chiral symmetry means that there exist a matrix which anticommutes with a 
Hamiltonian. 
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In this paper, we investigate in detail the phase diagram of the disordered fermion model with tt-Aux by means of 
the fermionic replica method to clarify these points. We show that this model, described by double Dirac fermions, 
has three kinds of chiral symmetry in special cases. Interestingly, these different chiral symmetries lead to different 
universality classes. To be concrete, an effective action of the Green functions has Sp(n)xSp(n) symmetry in one of 
them and U(2n) symmetry in the others, both of which are broken down to Sp(n). It turns out that corresponding 
Goldstone mode can be described by Sp(n)2 WZW model and U(n)/Sp(n) nonlinear sigma model, respectively. 

This paper is organized as follows. In the next section, we introduce a lattice model, which is the same as Fisher 
and Fradkin studied, and derive Dirac Hamiltonian in the continuum limit. Although the most generic model has 
only time-reversal invariance, symmetry is enhanced on three surfaces in the space of the coupling constants. We 
focus our attention to one of them with Sp(ri)xSp(n) symmetry (others have U(2n) symmetry) and give calcualtions 
in detail. Based on the continuum Hamiltonian, we derive an effective generating functional of ensemble-averaged 
Green functions in Sec. III. We show that Sp(n)xSp(n) symmetry is spontaneously broken down to Sp(n). In Sec. 
rV, we derive an effective action for the Goldstone mode, resulting in Sp(n)2 WZW model (U(2n)/Sp(n) nonlinear 
sigma model for the others). In Sec. 0, we present the phase diagram. In Sec. VI, we give a summary and concluding 
remarks. 



II. MODEL 



A. A lattice model and its continuum limit 



We define a baisc lattice model on the square lattice in two dimensions. For convenience, we introduce similar 
notations used by Ludwig et al ||26| in the following way. The square lattice in two dimensions is divided into unit 




FIG. 1. The square lattice on which the model is defined. Unit cells (plaquettes) are enclosed by the dotted-lines. Four sites 
on each plaquette are numbered counterclockwise. Thick lines means that nearest neighbor hopping is negative on them due 
to TT flux. 

cells, each composed of a plaqutte enclosing four sites. A plaquette is labled by a set of integers j = {jx,jy) and four 
sites on it are numbered 1,2,3 and 4 counterclockwise, as illustrated in Fig. |l| With such a notational convention, a 
pure and disorder Hamiltonian denoted respectively as Hp and iJj are defined by 



j a=l 



(2.1) 



where x = (1,0), y = (0, 1), and v'^j is random real number due to nonmagnetic impurities. Negative signs of the 
third and the sixth terms in Hp are due to tt flux imposed. Now let us construct an effective Hamiltonian near the 
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band center in the case of weak disorder. To this end, we first derive the continuum Hmit of the pure model, and next 
we take into account the effects of disorder. 

Fourier transformation Caj — J (^J^-i e'^^-'Ca{k) for a = 1, ...,4 leads to 



1 + e^^" 1 + e^^^y 



c{k), 



where c^(fc) = {c\{k), ...,c\{k)). It is readily seen that the Hamiltonian gives zero-energies at fc = (tTjTt). In the 
vicinity of them, the model can be effectively described by Dirac fermions. To see this, set k — (tt, tt) + aop with the 
lattice constant ao, expand the Hamiltonian to the first order of p, and introduce a field operator in the continuum 
limit CQ((7r,7r) + aop) = Og ^tpaip) for a = 1, ...,4. Then we have 

J^:^i^Hp) i-<^2 «> IPx - CTi (g) a2Py) i^{p), 

where V'^ = ("017 "02' V"-!; V'l)- Accordingly, effective Hamiltonian near the zero energy can be expressed by the following 
field theoretical Hamiltonian 

Hp = J d^xi;\x)Hp^{x), Hp=j^td^, (2.2) 

where x = aoj, ip{x) = J ^^J^i e^P^ip{p), and 7 matrices are defined as 71 = 172 (X" 1 and 72 = ci (g) tT2- 
By using the relation between the lattice operator Caj and the continuum field operator ipa{x), 

the continuum limit of the disorder potentials are readily derived; Hd — J d^x^^^^^ ^\(x)v'^{x)ilja{x) , where v'^{x) is 
the slowly- varying part of v'^j defined as v'^^j = aov'^{x). Taking suitable linear combination of v'^, we can rewrite the 
random potentials as 

r * 

Hd= d'^xij''{x)ndij{x), Hd = Y,vAx)a^. (2.3) 

i/=i 

where ai = 1 (g) 0-3, 02 = (T3 g) 0-3 , 03 = (T3 g) 1 and 04 — 1 (g 1. We assume that the four disorder potentials obey 
Gaussian probability distribution 

P[v^{x)](xe~^'''''''\ for i/ = l,...,4. (2.4) 
It may be more convenient to switch into another basis via suitable rotation for i/i, 

71 = (72 g) 1, 72 = 0-1 g) 0-3, 

ai = lg)cr2, a2 = CT3 (8) ^2 , a3 = 0-3 g) 1, 014 = 1 g) 1. (2.5) 

In what follows, we refer to the space where these matrices live as V. It is distinguished from the replica space 
Wr introduced momentarily. The total Hamiltonian we investigate in this paper is H = Hp + in Eqs. (2.2) 
and (2.3) with (2.5). It is noted that this is the same model as Fisher and Fradkin studied in the case 171 = 
92 = 93 — 34- They showed that the effective field theory is 0(2n, 2n)/0(2n)xO(2n) nonlinear sigma model for the 
bosonic replicas [or Sp(2n)/Sp(n)xSp(n) nonlinear sigma model for the fermionic ones], concluding thereby that all 
states are localized. According to Zirnbauer, this class is referred to as All Although a slight relaxation of the 
condition 91 = 92 ^ 93 ~ 94 still leads to the same universality class, several subspaces in {gi, 92, 93, gi) have not 
only time-reversal symmetry but also additional enhanced symmetry specified below, and therefore belong to different 
universality classes. 
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B. Symmetries of the model 



In this subsection, we examine symmetries of the present model, which plays a crucial role in the classification of 
the universality classes for disordered systems . First of all, the most general Hamiltonian has unbroken time-reversal 
symmetry. In the original basis, it is expressed by TC* = Ti, whereas in the rotated basis where the Hamiltonian is 
described by the matrices in Eq. ( ^.5[ ), it is expressed by 

n^Tn''T-\ r = 1(8(71. (2.6) 

If the model does not have any other symmetries, the universality class is ^11 and hence all states are localized. There 
are several subspaces in (gi, (72, ff3, 54), however, on which symmetry is enhanced. To be concrete, the Hamiltonian 
( ^.2[ ) has following symmetry 

T'l = 0-3 (g) 0-3, for 53 54 = 
n = -VjHV-\ with { V2^(Ji® (71, for 52 = .g4 = . (2.7) 

Vz^ CFi® z(72, for 51 = 54 = 



Recently, Hatsugai et.al |9[| proposed the same model as Eq. {'2A_) but with random hoppings instead of random 
on-site potential, which we will refer to as HWKM model. They have shown that such Hamiltonian has not only 
time-reversal symmetry but also another one often called chiral symmetry. It has been shown that this 

symmetry is responsible for the existence of delocalized critical point at the zero energy. Chiral symmetry in this 



case is expressed as Ti — —ViHV^ with 7^4 = as (8) 1, similarly to Eq. (|2.7D, and the universality class is specified as 



BDl or in other words U(2n)/Sp(n) (Gl(2rt,R)/0(2n)) nonlinear sigma model for fermionic (bosonic) replica method. 



Considering these, it is easily guessed that the symmetry (2.7) also plays a crucial role in Anderson localization for 



the present model. In what follows, we refer to the symmetry in Eq. (2.7) also as chiral symmetry. 



III. REPLICA METHOD 

It is so tedious to repeat similar calculations that we concentrate on the case Vi, and hence we often denote Vi 
just as V for short. But with slight modification, any Vj can apply to the formulation. At the moment, the difference 
between Vi and the others is whether Vj commutes or anticommutes with T and whether it is real symmetric or 
antisymmetric. Surprisingly, this difference, which seems trivial at first glance, leads to a different universality class. 

In order to study the trasport properties of the model near the zero energy, we need to compute the Green functions 

G[±ie) ^ 



±i€-n 



With the V symmetry, however, we do not need to distinguish the retarded and advanced Green functions, since they 
are related with each other via the relation |^ , 

G(+ie) = V^—V-^ = -VG{-ie)V-\ 
le+rL 

The Green function G{ie) can be computed by a replicated functional 



Z ^ 



with 



£ = ^ ^ 'ip*^{ieli - 'H)jktpka = trM/^'0^(«el4 - 'H)^', (3.1) 

a=l j,k=l 

where the subscript a denotes the n species in the replica space Wy, ■0^^- = tpjai ^^'^ the usual convention of the matrix 
product is applied to the fermion matrix field ip*^ and ipja- Although such a notational convention is not always 
necessary, it is quite convenient for us to save many complicated subscripts. It should be noted that the fields tAqj 
and Tpaj are completely independent Grassmann variables. Therefore, they obey in general different transformation 
laws. 
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A. Auxiliary space 



Based on the Lagrangian (3.1), we will construct an effective theory of order parameter fields for the replicas by 
averaging-over disorder. In this process, it is quite important to fully incorporate the symmetries (2.6) and ( ^.7] ) in 
the space V. To this end, we introduce an auxiliary space to reflect these into the replicas, according to the method 
developed by Zirnbauer Namely, we extend the replica space as ^ W = Wr Wt Wc, where Wt (Wc) 
stands for the space associtaed with time-reversal (chiral) symmetry. We then introduce a 4n x 4 (4 x 4n) matrix 
field (^'ja), 



* = (x+,x-); x± = x^T' x±^'P±x, 



where 



1 

71 



V2 



(3.2) 



(3.3) 



V± = ^3^, and -0' stands for the transpose of ■(/>. Since the components of these fields are not independent of each 
other, the fields are subject to the constraints 



$ = -TT^-P, * = ■P^'TT, TT = 1„ 



1(72 I2, 

1 12 (8) 0-3. 



(3.4) 



Note that the matrices 7 and tt are defined in W, whereas T and V are defined in V. We immediately find 

=T(*$)*r"\ = -7?(^'$)7?-i. (3.5) 

These relations tell that the matrix '5^ in V has the same symmetries as TC has. It is also easy to find alternative 
relations for the matrix \E'\E' in W, 



*\I' = 7(vl/\E')'7 \ vj/v]/ = -7r(vI/\I/)7r- 



(3.6) 



These relations actually specify the symmetry of the model, as we shall see lator. 

Although the potentials ai and a2 have time-reversal and chiral symmetries, the converse is not true: The symme- 
tries (2.6) and (p.7D admit other potentials /3i = (Ti ® 1, /?2 = o'2 (8) 0-3, /33 = 1 ® cti and /34 = (T3 ® cri. The reason why 



there are no beta matrices in Eq. (2.3) is that the basic lattice model (2.1) includes only on-site random potential. 
Namely, if we take all kinds of possible disorder for the lattice Hamiltonian including random (nearest neighbor and 
diagonal) hoppings as well as random on-site potentials, we h ave a continuum Hamiltonian with all kinds of possible 
disorder potentials specified by the relations (2_^) and (2/7). Such a model may be referred to as a model with 
maximum entropy. Since we are interested in the model with on-site random potentials, we have to remove these /3 
matrices. For this purpose, we take account of additional conditions 



Hd = -THdT, Hd = SHdS, f = (73 ® 1. 



(3.7) 



It is readily seen that the potentials including beta matrices do not satisfy these relations. Let us then introduce the 
fields 



(3.8) 



We refer to this extra space as W° . It is easy to verfy that the matrix \I>°\I'° obeys the same relations that Ti^ does in 
(3.7) and in (|3.5| ) with extension 7 ^ 7 (8) I2 ® and tt ^ tt (g) I2 ® (T3, whi ch m eans that the undesirable potentials 
are in fact excluded. By using these fields we can reexpress the Lagrangian (3T) as 



where a; = 1, 



I2 ® <J\ and ^ = 7^5^,. 



(3.9) 
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B. Effective Lagrangian averaged-over disorder 



Now we are ready to ensemble-average the Lagrangian ( |3.9D over disorder. In this subsection, we assume 171 = .92 = 5 
for simphcity. More general case 171 7^ 52 is briefly discussed in Sec. [II D| . The probability distribution (2.4) can be 
expressed as 



(3.10) 



Therefore, the Lagrangian is converted into 



1 



2.9 

= —iYv (Ud + g-^°^°) - -trv(*°$°)^ + trvF(iet^$* - ^i^*), 
2g \ / 2 

where we have used a relation tj:w(g>w°^°T^d'^° = — trv'7Yd^'°^'°- It should be noted that integration over TYd is 
readily performed, since the matrix vI'°v]/° has completely the same symmetries as Hd has. This is the reason why we 
have introduced several auxiliary spaces in the last subsection. By using try(^'°5'°)^ — — trvi/®vK° (^°^°)^ again and 
the definition of (B.8), the gaussian integration over Hd leads to 



C = tr^v{ieuj^^ - ^i^-^) + |- ^ trvy(5'C'^4')2, 



(3.11) 



T = l 



- We mention in passing that if the model have maximun 



where d = 1,02 = 5, ©3 = i^, O4 = iT£, and g' ^ 
entropy, there appears only one interaction term with Oi = 1 and g' = g. Eq. ( [3.11 ) is an effective Lagrangian which 
generates the emsanble-averaged Green functions for the present model. It is easy to confirm that new matrices are 
subject to 



^Or"^ = j{'90r'^Yj^^, *Cr* = T7r(*C'r*)7r"^ for 



T=l,2 

T = 3,4, 



(3.12) 



So far we have obtained an effective action including four fermi interactions by averaging over disorder. The next 
task is to introduce auxiliary (order-parameter) fields to decouple them. To this end, we add the following terms of 
auxihary fields Qr (t = 1, 4) to the Lagrangian ( [3.11 ): 



1 4 

= - ^ ^ tr,v(Qr + g'Wr^ + Srlteujf 
■J 1 — 1 



Then, the Lagrangian becomes 



(3.13) 



(3.14) 



Eq. (3.13) requires that Qr should obey the same symmetrirs that ^O^-^' does in Eq. (3.12). Therefore, 



for 



T=l,2 

T = 3,4. 



(3.15) 



The antihermiticity ensures that the integration over Qr in Eq. (3.14) converges. 

This Lgrangian (3.14) has a global Sp(n)xSp(n) symmetry. To see this, let us onsider a transformation 5* — > '^g^^, 
'i' g'i/, and Qr gQrg^^- Since this transformation should keep the relations (3.15), g is subject to g^'g — 
1:5*75 = 7 ^-nd g~^TTg = tt- Therefore, g has the following form 



5+ 



g- 



, 5±5± = 1, 5±7o5±=7o, 7o = ln«'icr2, 



where 2n x 2n matrices g± belongs to Wrt = W^r ^ Wt and the explicit matrix above represent the space Wc- The 
conditions for g± require g± € Sp(n). 
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Next let us solve Eq. ( |3.15D , 



Qr= (^_^^t ^'^J' for r=l,2, 
Qr- (^^^-2 A- 2)' ^ = 3,4, 

where 2n x 2n m atric es Mr and are subject to Mr = —joMr^^^^, = 70^^*7^^'"'^, and A^^ — —A^. Using the 
definition in Eq. (3_^), the Lagrangian is now written as 

1 ^ i 

2 

- ^ trvy^, (M,X-OrX+ - Mlx+OrX- + A+x+0,+2X+ + • (3.16) 

T = l 

The tranformation laws under g are 

X± ^ XisL X±^5±X±, Mr ^ g+Mrgl, ^ giA^g^. (3.17) 

So far we have examined the case Vi. It has been shown that the Lagrangian has Sp(ri)xSp(n) symmetry and hence 
the model belongs to CI, which is the same class as d-wave superconductors with unbroken time-reversal and spin- 
rotation symmetries. Contrary to this, in the other cases associated with V2 and V3, symmetry of the Lagrangian is 
quite different. Similar calculatio ns ab ove show that these Lagrangians have U(2n) symmetry insted of Sp(ri)xSp(n) 
dut to the fact that 7 defined in {^A) is 7 = 1„ (8) ia2 ® (Ti if we take the same tt. This is the same class as HKWM 
model, which is often referred to as BDl or chiral orthogonal class. 



C. Effective potential 

Based on the Lagrangian in the last subsection, we derive the effective potential to study the spontaneous symmetry 
breaking. We assume that the order-paramters are diagonal in the space Wi-. Namely, set 

Mr = i diag(mT-i, ...,to^„) ® I2, Af = i diag(a^i ± bri, ■■■,arn ± bm) t = 1,2, 



where mra{x), araix) and bra{x) are real fields. Then the Lagrangian ( 3.16 ) becomes 



2 2 " 2e 



^ r=l Q = l ^ Q=l a=l 



where is 4 x 4 matrix defined by 

Va = -~i{miai §51 + m2aCrz ® 1) + OiqI ® (Ti — ibi^a^ (g) (T2 + a2aCr3 (g) (Ji — ibi^l (g) (72. 

Integrating out the fermi fields, we have an effective potential Va{{m, a,b}) — S{{m{x) — m,a{x) — a,b{x) — 
&})/volume of one- loop order, 

2 ^ ^ C cPk 2e ^ 

Va({m,a,6}) = — X! ("^ra + aTQ + bla) - / -y^lndety(/^-h Wa) + — X! "^i"- 

Since the potential Va is a function of six variables, it is hard to find out the true minimum. However, it is readily 
expected from the symmetry breaking term that toiq has a finite value. Assuming it, we have tried to find out 
nontrivial sulution dVa/drrira = dVa/dara = dVa/dbra = 0, but in vain. Although possibility of some other 
solutions still remains, the solution mia 7^ and others=0 is of course a natural solution. In this case, the saddle 
point equaiton for mia becomes 

^ _ dVg _ 4mia _ 2 ^ 2mia 



dniia g' J (27r)2 fc2 + 
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47r 1 _ 2-77 ~ 

which gives nontrivial solution mi„ = m = A(e"?" — 1)2 '-^ Ae This result suggests that the operator (^Ci^) 
is relevant whereas others are irrelaevant around the present saddle-point. Therefore, in the following analysis, we 
neglect the irrelavant interactions and take only the relevant interaction into account. In the next section, we will 
discuss slow fluctuation aroung the saddle point. But before it, we will see how the anisotropy gi 7^ (72 modifies the 
theory so far obtained. 



D. Anisotropy effects 

In the case where gi =52, we could treat them simultaneously in Eq. ( B-lOj ). Contrary to this simple case, we 
have to treat the two potentials ai and a2 separately in anisotropic case. To this end, first introduce a new matrix 
= CTi (g) I2 under which they transform 



We next introduce 



ai = Ta\T ^, 0.2 — —J-a2F ^ . 





Let us denote Tidj = '^jOtj for j — 1,2 and the extra auxiliary space as W. Since the matrix ^'j^'j {j — 1,2) has 
completely the same symmetry property as Tidj , we can average over disorder in the following way: 

^trvnli + ^tYvnl2-tyW(SWo^W' f$?Hdl*? + $2Wd2«'2 

^trvy,g,iv°®w' + ^trvv^vF°(8iw(^i*i)^ 

91+92, /',T,0,TrO\2 I 91" 92, /-.TfO -7-,TrON2 

= — ^ — trM/^vF°(* * ) + — ^ — trvi/®w=(* ) , 

where arrow menas that the integration over Tidj is carried out. This equaiton shows that anisotropy gi ^ 32 yields 
an additional term including J-. However, it turns out that this term is also irrelevant and therefore can be neglected, 
following similar discussions to the previous subsection. 



IV. EFFECTIVE ACTION FOR THE GOLDSTONE MODE 



The saddle-point Mi = 77il„ (8) I2 is invariant under 17+ = g- type transformation in Eq. (3.17) but not invarinat 
under g+ = g^_. This means that Sp(n)xSp(n) symmetry is broken to Sp(n). Therefore, fluctuation around the 
saddle-point can be incorporated by parametrizing the order-parameter field Mi as 

Ml — ^ZTO^ = iTO^^ — imll, 

where ^ and U G Sp(n), i.e., U^jqU = 70, U^U ~ l2n- The field U descirbes the Goldstone mode around the saddle 
point. Then, the Lagrangian becomes 

^ = -trM/« [x+i^X+ + X-i^X- + im{Ux-X+ + U^X+X-)] + ^t^w^^U + U^) + 

The last constant term will be neglected, since it vanishes in the replica limit n — ^ 0. In the following subsection, we 
construct an effective theory for small fluctuation U around the saddle point. 



A. Integration over fermi fields 

So far we have treated x ^^'^ X ^ matrix fields to save complicated subscripts. It has actually worked well 
especially for the discussions of the symmetry. However, in order to extract an effective Lagrangian for U via chiral 
gauge transformation introduced below, we cannot treat the space V and W separately any longer. Therefore, we 
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switch into the standard notation for fermi fields as follows: Note the following identity, for example, tTWrtUx-X+ = 

U^X+j where the last arrow means that we have made a 70 rotation to x and x and we have regarded x± as a vector 
with 8n components in 1/ (g) Wrt • Then the fermion part of the Lagrangian can be denoted as 

- = X+i^® l2«X+ +X-i^®'^2nX- (X-I4 U^X+ + X+U «) Ux-) 

= X ® l2« + {V+ +V-® U)] X- (4.1) 

Although this Lagrangian seems to be standard, the integration over x is not straightforward. It is because the 
elementary fields 1/'^ and 1/' are doubly included in x and x defined in Eq. (3.3). It is possible, however, to rewrite 
this Lagrangian as 



i{d2 +«ai)i2„ 



i{d2 ~ idi)l2n 
imU 



(4.2) 



where ?/;' and ip' have An components whose explicit subscripts in the above definition are those for V . It should be 
noted that this Lagrangian is described by independent field variables ip' and ijj' only. As a result, we can integrate 
out fermi fields, which results in 



= Dety^vKr 
= Detf 



imW i{d2 ~ idi)l2n 
i{d2 + idi)l2n imU 



{r+ + V-® U)] . 

Based on this result, we derive an effective action for the fluctuation U in the next subsection. 



(4.3) 



B. Chiral gauge transformation 

Suppose that the fluctuation U{x) around the saddle-point varies slowly with x. Then, we can expand, with respect 
to derivatives, the determinant derived in the last subsection. To this end, make the folowing transformation 

x'i^U^X+, 

V — V V — \^ 

A — A — 1 A — A — • 

Here and hereafter, we often denote I4 (8) J7 just as U for simplicity. Without confusion, we also use e.g., V instead of 
'P <8> l2n. The Lagrangian is now converted into 

- Cf ^ X^ {i P + im)x^ , 

where 

D^ = D+f, + d-^; D+^ = D^V+, = (4.5) 

with 

Integration over fermi fields xF and x^ gives also a determinant, and the action for U is 

S^-\ InDetv^w,, {ip + im) - In ./([/) + ^ ^ d^xtrw., (U + U^) (4.6) 



where J{U) is a Jacobian due to the transformation (4.4). It turns out that expansion of the determinant yields 
no second order terms with respect to the derivative, and hence we neglect the contribution from it in the leading 
order approximation. In what follows, we then concnetrate on the calculation of the Jacobian. We mention in passing 
that we will use the Lagrangian (4.1) to compute the Jacobian, although (4.2) is less confusing. 
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Let us consider the following eigenvalue problem for the Dirac opeartor (4.5) to define and regularize the Jacobian 
of the path integration, 

where right- and left-eigenvalue equations are introduced, since the Dirac operator is non-hermitian. The eigen- 
functions are normalized via the inner product {<j)l^,(pm) = J dl^xcfx^Lpm = Snm- Assuming the completeness of the 
eigenfunctions, the fields x and x are expanded as 

The measure is then defined as 

Vx"" ^HdaH or I?x^=J]d6^. 

n n 

It should be noted that if one uses the Lagrangian (^^) described by independent field variables ip' and ip' only, 
the measure is defined by Tyip'Vip' as usual. In the present case, however, either Vx or Vx is enough, since each 
already includes all independent variables. It seems to be unusual, to be sure, but we follow the calculation with the 
latter notation. Using the definitions for the path integration measure, we can compute the Jacobian for infinitesimal 



small change for the gauge transformation (4^), x^ ~* t^u+su _ _j_ ^^u ^^^^ ^u+su _ _^ g^u ^j^^jj 



Sx'^ = -r+ ® WdUx'^ and dx'^ = x'^V^ ® WSU as follows 

m 

where 

Nmn = -{(blr-®UUU,ipn). 

Therefore, the measure changes as '^"n = ^ety^w^t^Mi^ — -^)n^"n^'^^ and H^^n = ^ety^Wrt^Mi^ — 
TV) ]^ , where M stands for the space for the eigenmodes of the Dirac operator. Noting that det(l -I- eX) ~ 

gctrX+oie ) ya,ii(j for small e, we reach 

Vx"" ^Vx^'+'^'e''^-, 

where 6T^ = Tiy^w^t^MN = J (fxY^n^'^v^WrtNnn and i5r~ = Trv^w^t^AiN = J d^xX]„ try^w-.t A'^nn. Here TrM 



stands for the trace for x and M. In principle, we can choose either or T~ as the Jacobian in Eq. (4^) if they can 
be regularized. In what follows, we will apply standard regularization scheme in Refs. [^,^. Nevertheless, it turns 
out that they still give rise to divergent terms for each F. Fortunately, as we shall see below, they cancel each other 



in Sr^ + ST~. Therefore, it is suitable to define the Jacobian in Eq. (4.4) as a square root of the product of the naive 

1 ~ 
Jacobians as Dety^^ ®m(1 — N){1 — N). Namely, 



s\nJiu)^sr = -{sr^ + ST~). 



It is possible to show that the Jacobian based on Vip'Tyip' for the Lagrangia n ([4.2 ) indeed yields the same Jacobian. 



We now give the outline of how to calculate the Jacobian following Refs. [^8[|2S|]. Using the definition of the Tym, 
we have 
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n 



where ultraviolet cut-off A is introduced, which will be set A ^ c» at the end of the calculation. Using the relation 

= p^^_ + ^_p_^_, we have 



1 /■ d^k 



2 y (^Tr.,.^^^-"^t^^^t^ 



ikx 



where in the last hne we have switched the basis for M from the eigenmodes of the Dirac operator into plane waves. 
Expanding A and taking the limit A ^ cxj, we reach 

= i y" -0^e-''"TTv»w., [V ® U^5U{id xD) + l® U^5U{d ■ D)] . (4.7) 
where a x b = fT^j/fl^bi^ and a ■ b — with 

CTmi' = ^bn'Ti'] = £^'''^3 «) era, g^u ^ = Sf,i,l(»l. (4.8) 
Therefore, taking trace in x- and 1^-space, we have 

It is nowadays well-known that the action function satisfying above includes not only the usual kinetic term but also 
the WZW term, 

r = 5{Vzw = ^/ d^xtrw^^d^Ud^U^ ^tLJ d^^e^->^ti'w^^U^d^UU^d,UU^dxU, (4.9) 

with k = 2. Adding the symmetry-breaking term in ( [4.6| ), we end up with an action 

S = + ^ / rf'^^trvi.,. ([/ + U^). (4.10) 

Therefore, it turns ou t th at the Goldstone mode of the model wi th ch iral symmetry of type Vi can be described by 

Sp(n)2 WZW model (O) with a symmetry breaking term as in ( 4.10 ). 

So far we have considered the case where Hamiltonian has the chiral symmetry specified by Vi in Eq. (2.7). As to 
the other cases, similar calculations lead to U(2n)/Sp(n) nonlinear sigma model, which is so-called chiral orthogonal 
class or BDl, the same universality class as HWKM model. 



V. PHASE DIAGRAM 



We present a schematic phase diagram at the zero energy in Fig. |[ which is valid for weak disorder, i.e, for small 
gj. In the case gs = g^ — 0, which is denoted by a white surface in the figure, the model belongs to class CI and is 
described by Sp(n)2 WZW model with some irrelevant perturbations discussed in Sec. [II C and HID. The reason 
why the present WZW model is level-2 is that the continuum limit (2.2) leads to double Dirac fermions. One of the 
main origin of the irrelevant interactions is that the lattice model has only on-site potentials without any hopping 
randomness. Namely, the model does not have maximum entropy. Another origin is the anisotropy gi 7^ (72, as 
discussed in Sec. HID. Besides, the continuum model should have basically many higher order irrelevant interactions. 



Therefore, at the zero energy any points on this surface flow to the Sp(n)2 WZW fixed point, denoted by a black 
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g 




gl=g2=g3=g4 



FIG. 2. Schematic phase diagram at the zero energy in the case of weak disorder. We omitted the axis 34. Black circle 
indicates the Sp(n)2 WZW fixed point. Shaded surfaces stand for "critical surfaces" described by U(2n)/Sp(n) nonlinear sigma 
model. Any points away from these three surfaces flow to strong coupling limit of 0(2n, 2n)/0(2n)xO(2n) nonlinear sigma 
model. 



circle in the figure. It turns out that around the fixed point the density of states obeys the scaling law p{E) ^ Ett ^ 
following from the exact results obtained by Knizhnik and Zamolodchikov . 

On the other hand, the universality class of the two surfaces (71 = 54 = and (72 = ff4 = 0, shaded in the figure, is 
different. They belong to the class of U(2n)/Sp(n) nonlinear sigma model or in other words, chiral orthogonal class 
or class BDl, which is the same universality class as that of HWKM model. It is well-known that the beta function 
of the renormalization group vanishes and the coupling constant of the nonlinear sigma model is not renormalized. 
Hence any points on these surfaces are critical. The density of states is conjectured to be divergent at the zero energy, 
according to p{E) ^ ^ exp[— a-\/ln(A/i?)] with a positive constant a 

These surfaces, one CI and two BDVs, are unstable: Out the surfaces renormalization group fiows are toward the 
strong-coupling limit of Sp(2n)/Sp(7i)xSp(n) nonlinear sigma model, as shown by Fisher and Fradkin for gi = g2 = 
93 = 54 case. Therefore, all states away from the surfaces are expected to be localized. 



VI. SUMMARY AND CONCLUDING REMARKS 



We have studied the phase diagram of the disordered fermion model with random on-site potential defined on a 
two-dimensional square lattice with 7r-flux. The continuum limit near the zero energy yields Dirac fermions with 
random potentials specified by four independent coupling constants. It has been argued that in addition to time- 
reversal symmetry, this model has so-called chiral symmetry in three subspaces in the space of the coupling constants. 
Here, chiral symmetry means that there exist a matrix anticommuting with the Hamiltonian. What is interesting 
is that this model alows three kinds of chiral symmetry, which is due to the species doubling of the lattice fermion. 
It has been shown that one of them leads to Sp(n)xSp(n) symmety of the Lagrangian for generating functional of 
emsemble-averaged Green functions whereas the others lead to U(2n) symmetry, both of which is broken to Sp(n). 
It has been found that the fluctuation around the saddle point is described, respectively, by Sp(n)2 WZW model for 
the former and U(2n)/Sp(n) nonlinear sigma model for the latter. Based on these results, we have proposed a phase 
diagram. 

The Dirac fermion model with Vi symmetry calculated in detail in the text is quite similar to that for the disordered 
d-wave superconductors with spin-ratation and time-reversal symmetries. Actually, effective Lagrangians for these 
systems have Sp(n)xSp(n) symmetry and hence both the models belong to the same universality class CI. However, 
fixed-points are quite different due to the WZW term. It was shown jl^,^ that in the d-wave superconductors the 
WZW terms for several species of Dirac fermions cancel each other and the resultant theory predicts localization of 
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quasiparticles at the band center. Contrary to this, the model treated in this paper is described by level-2 WZW 
model, which means that it has delocahzed states. The difference lies in Eq. (4.8): In the present model the WZW 
term survives owing to the relation (t^j, = [7^,7i,]/2i — e^yVx, which is due to V\ = 17172- This is just the same 
definition of so-called 75 in the particle physics literature, although the present 7 matrices are not 2x2 but 4x4. 
Conversely, we can say that chiral symmetry only in this sense yields the chiral anomaly and hence the WZW term, 
although we often use the same words for other similar symmetry in the Anderson localization literature. Considering 
the fact that the lattice model having Vi chiral symmetry is realized in the case v'^ = —v[ and v'^ = —v'2 in Eq. (^.l|), 
which is indeed a fine-tuned model, we conclude that the delocalization due to the existence of the WZW term seldom 
occures in nature. 
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